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Abstract 

With perfectly balanced gain and loss, dynamical systems with indefinite damping can obey the 
exact "PT-symmetry being marginally stable with a pure imaginary spectrum. At an exceptional 
point where the symmetry is spontaneously broken, the stability is lost via passing through a non- 
semisimple 1 : 1 resonance. In the parameter space of a general dissipative system, marginally 
stable T^T-symmetric ones occupy singularities on the boundary of the asymptotic stability. To 
observe how the singular surface governs dissipation-induced destabilization of the "PT-symmetric 
system when gain and loss are not matched, an extension of recent experiments with "PT-symmetric 
LRC circuits is proposed. 

PACS numbers: ll.30.Er, 03.65.-w, 41.20.-q, 45.20.-d, 46.40.Ff, 45.10.-b, 02.40.Vh, 45.30.-t-s, 47.20.-k 
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Introduction. The notion of PT-symmetry entered modern physics mainly from the side 
of quantum mechanics. Parametric famihes of non-Hermitian Hamiltonians having both 
parity {V) and time-reversal (T) symmetry, possess pure real spectrum in some regions of the 
parameter space, which questions need for the Hermiticity axiom in quantum theory [l-4|. 
First experimental evidence of "PT-symmetry and its violation came, however, from classical 
optics in media with inhomogeneous in space gain and damping and electrodynamics 

PT-symmetric equations of two coupled ideal LRC circuits, one with gain and another 
with loss, have the form 

z + Dz + Kz = 0, (1) 

where dot stands for time differentiation and the real matrix of potential forces is K = > 
while the real matrix D = of the damping forces is indefinite Q]- 
For the problem considered in 3], we assume that 



D = D 



VT 




K = K 



VT 




(2) 



= {zi,Z2), and 6, k and k are non-negative parameters. Eigenvalues of TivT have equal 
absolute values and differ by sign, indicating perfect gain/loss balance in system ([T]) with 
matrices ([2]). The coordinate change Xi = Zi + iz2, X2 = xl, x^, = Xi, and x^ = X2, where 
i = and the asterisk denotes complex conjugation, reduces this system to «x = Hx, 

where the Hamiltonian 

/ i 0\ 

i 
—ik K, i6 



H 



(3) 



y —K —ik i6 I 
diag (1,-1, -1,1)) 



is PT-symmetric (PH*=HP, P= 

In real electrical networks, additional losses may result in the indefinite damping matri- 
ces that possess both positive and negative eigenvalues with non-equal absolute values. A 
systematic study of dynamical systems ([1]) with such a general indefinite damping, has been 



initiated in 
biology 



10 



12|-Qj. In 



nj in t 



15- 



le context of distributed parameter control theory and population 



17j gyroscopic stabilization of system ([T]) was considered, because 



negative damping produced by the falling dependence of the friction coefficient on the sliding 



2 



velocity, feeds vibrations in rotating elastic continua in frictional contact, e.g. in the singing 



wine glass 



18l-l2l| . In 22| a gyroscopic PT-symmetric system with indefinite damping was 



shown to originate in the studies of modulational instability of a traveling wave solution of 
the nonlinear Schrodinger equation (NLS) j^l- In nonlinear optics, a challenging problem 
of stability of localized solutions (solitons) is related to the indefinite damping, because 



stable pulses in dual-core systems frequently exist far 



perfect matching of gain and loss (PT-symmetry) 



rom the conditions that provide a 



25| . Recent techniques proposed 



for the stabilization of the solitons in two coupled perturbed NLSs include introduction of 



PT-symmetric nonlinear gain and loss [26i] which signs can be periodically switched |27l. |28|. 
Therefore, indefinite damping is a basic model to study how a localized supply of energy 
modifies the dissipative structure of a system fl^ . 

In general, the eigenvalues (A) of system ([1]), when it is assumed that z ~ exp(At), are 
complex with positive or negative real parts corresponding either to growing or decaying in 
time solutions, respectively. Asymptotic stability means decay of all modes. 

A two-dimensional system ([T]) with D = 61) is asymptotically stable if and only if trD > 

(trKD - (Ti(K)trD)(trKD - (72(K)trD) 



and < 5^ < (5^^, 



5^ 

cr 



det DtrDjtr 



KD - trKtrDl 



(4) 



where cri(K) and cr2(K) are eigenvalues of K 10|,l29|. However, when simultaneously trD = 
and trKD = 0, the spectrum of the system ([1]) is Hamiltonian, i.e. its eigenvalues are 
symmetric with respect to the imaginary axis of the complex plane HI]. They are pure 
imaginary and simple {marginal stability) if and only if 6"^ < 6p-j-, 

1/2 



S-pT 



ai(K)- Va2(K) 



-detD 



(5) 



How the marginal stability domain of a indefinitely damped PT-symmetric system relates 
to the domain of asymptotic stability of a nearby dissipative system without this symmetry? 
The answer is counterintuitive already for the thresholds (jl]) and ([5]). Our Letter describes 
mutual location of the two sets, thus linking the fundamental concepts of modern physics: 
PT-symmetry jl-4| and dissipation-induced instabilities 30-32|. 

A potential system with indefinite damping. First, we extend the model ([T]) with matrices 
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FIG. 1: (a) In the half-space X > of the {ki,X, Y) space, where X = 61+62 and Y = 61—62, a part 
of the singular surface locally equivalent to the Pliicker conoid of degree n = 1, bounds the domain 
of asymptotic stability of system ([T]) with matrices ^ and n = 0.4 and /c2 = 1; PT-symmetric 
marginally stable systems occupy the red interval of self-intersection with two exceptional points 
(EPs) (black dots) at its ends, (b) The top view of the surface. 



([2]) by choosing the matrices of damping and potential forces in the form 



D 





^ 




( ki 












S2 1 






k2 / 



(6) 



where parameters can take arbitrary positive and negative values. For asymptotic stability 
it is necessary that trD > and det K > 291]. 

Introducing the parameters X = 5i + 82 and F = 5i — ^2, we use the Routh-Hurwitz 
stability threshold (jl]) where one should equate the right hand side to unity and replace the 
matrix D with that given in Eq. ([6]). The result is a quadratic equation for ki. Expanding 
ki{X) in the vicinity of X = 0, yields a linear approximation to the threshold of asymptotic 
stability in the (fci,X) plane 



ki=k2+^y 



y2 ^ 



VT 



V2 V+ ^ 



(7) 



Y^j- = 2 ^a/o"2(K) ± a/o'i(K) j , where ai = k2 — n and a2 = k2 + n are eigenvalues of 
the matrix K from Eq. in which ki = k2 that happens when X = 0, i.e. 61 = —62- 
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FIG. 2: Evolution of eigenvalues of system ([T]) with matrices ([6|) where k = 0.4 and k2 = 1. (a) 
The loops of pure imaginary eigenvalues (dark grey) between the EPs marked by green dots imply 
marginal stability of the "PT-symmetric system corresponding to A;i = 1 and X = 0. (b) Unfolding 
the EPs of the unbalanced dissipative system with ki = 1.2 and X = 0.2 and (c,d) its growth rates 
as functions of Y (red curves). The growth rates vanish at the lower values of Y not converging 
to the locations of the EPs of the "PT-symmetric system (red dots on a green curve) when /ci — )■ 1 



and X — >• along a ray in the {ki,X) plane (the destabilization paradox 



2g 



)• 



Therefore, on the line defined by the equations ki = k2 and X = in the {ki,X,Y) space, 
system ([1]) with the matrices ([6]) is reduced to the PT-symmetric system with matrices ([2]) 
that is marginally stable on the interval —Y^j- <Y< Y^^, cf. Eq. 

In Fig. [TJ^a) the vertical red line denotes this interval with Y-j^j- 0.817 calculated for 
^2 = 1 and K = 0.4. Along it PT-symmetry is exact, i.e. eigenvectors are also VT- 
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symmetric [l-4|. Hence, the spectrum is pure imaginary, see Fig. |2l The ends of the 
interval are exceptional points (EPs) [33] corresponding to the merging of a pair of pure 
imaginary eigenvalues into a double one with the Jordan block. Passing through these 
points of the non-semisimple 1 : 1 resonance with the increase of \Y\ is accompanied by 
the spontaneous breaking of the "PT-symmetry of eigenvectors although the system still 
obeys the symmetry. This causes bifurcation of the double pure imaginary eigenvalues into 
complex ones with negative and positive real parts and oscillatory instability or flutter when 
Yj^r^ < < Y^j 1 see Fig. [2]^a). The bifurcation at = Y^j makes all the eigenvalues 
real of both signs (static instability or divergence). 

What happens with the stability near the red line in Fig. dJ^a)? Fig. MJo) shows that, e.g. 
at the fixed ki = 1.2 and X = 0.2, the eigencurves connected at the EPs with Y = ^Y^j- in 
Fig. |2]^a), unfold into two non-intersecting loops in the (ReA, ImA, Y) space, manifesting an 
imperfect merging of modes 3J] owing to gain/loss imbalance. 

Now the stability is lost not via the passing through the non-semisimple 1 : 1 resonance 
but because of migration of a pair of simple complex-conjugate eigenvalues from the left- 
to right-hand side of the complex plane at \Y\ < Y^j- ~ 0.817. For example, tending 
the parameters to the point (1,0) in {ki,X) plane along a ray, specified by the equation 
X = fci — 1, we find that the thresholds of asymptotic stability converge to the limiting 
values of Y+ ~ 0.615 < 0.817 and Y_ ~ -0.531 > -0.817, see Fig. Wic4). The limits vary 
with the change of the slope of the ray. Therefore, infinitesimal imperfections in the loss/gain 
balance and in the potential, destroying the PT-symmetry, can significantly decrease the 
interval of asymptotic stability with respect to the marginal stability interval. 

Such a paradoxical finite jump in the instability threshold caused by a tiny variation in the 
damping distribution, typically occurs in dissipatively perturbed autonomous Hamiltonian 
or reversib 



namics 



36 



e systems 



29 



35| of structural and contact mechanics 



30 



32 



34| and hydrody- 



38[, as well as in periodic non-autonomous ones [39|. We have just described a 



similar effect when the marginally stable system is dissipative but obeys PT-symmetry. 

A reason for the dependence of the limiting critical value of Y on the direction of approach 
follows from the linear approximation ([7]), which defines two straight lines orthogonal to the 
y-axis. When Y changes from —Y^j- to Y^^, the straight lines ([7]) rotate around the Y- 
axis. We remind that a set of points swept by a moving straight line is called a ruled 
surface 40|, |4l|. A right conoid is a ruled surface generated by a family of straight lines 
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that all intersect orthogonally a fixed straight line (the Y-axis in our case). Therefore, 
Eq. ([7]) defines a right conoid in the (fci, X, F)-space. In order to identify its type, we 
observe that Eq. ([7]) results in a cubic equation for Y. The third-degree term in it can be 
neglected when \Y\ < Y^^. Resolving the remaining quadratic equation and introducing 
the polar coordinates (p, 0) in the {ki, X) plane as ki = k2 + p cos and X = we find 



a parametric surface 



(p, 0) ( ^2 + PCOS0, 



psm( 



2k 



: Sm I 



(8) 



'k^ Vk2 

This is a canonical equation for the special type of the right conoid known as the Pliicker 
conoid of degree 1 — a singular surface with one horizontal and one vertical interva 



intersection 



40 



41| . The latter has at its ends two Whitney umbrella singularities 



42 



se. 



44|. 



f- 



Near the interval — Yj^j- <Y< Yj^j- shown in red in Fig.[T]^a), the boundary of asymptotic 
stability given by Eq. (j4j) converges to the Pliicker conoid (JTj), which is its exact linear 
approximation. The latter, in turn, is approximated by the ruled surface ([H]) that is in a 
canonical form for the Pliicker conoid. Qualitatively, all the three surfaces have the same 
singularities visible in Fig. [1] 

The approximation of type (|8]) can also be obtained from the perturbation formulas for 
splitting double semi-simple eigenvalues ±ik2 (diabolical points) corresponding to k = 0, 



ki = k2 and 612 = 0, see 



20 



21| . The Pliicker conoid of degree 1 singularity on the 



boundary of the asymptotic stability domain generically occurs as a result of the unfolding 
of the semi-simple 1 : 1-resonance 29. 40|. 

The "PT-symmetric marginally stable system studied in [7|, occupies a common 'handle' 
of the two Whitney umbrellas on the Pliicker conoid surface. The surface forms an instability 
threshold for the nearby systems with the gain/loss mismatch and additional coupling in the 
matrix of potential forces. These imperfections are realizable in the physical LRC-circuits. 
This opens a way for the experimental investigation of dissipation-induced instabilities and 



related paradoxes that are common for very different dynamical systems [30|, |42|. Indeed, 
since the singular geometry behind the destabilization paradox in dissipatively perturbed 
Hamiltonian, reversible, and PT-symmetric systems is the same, the experiments with the 
near-PT-symmetric LRC-circuits promise to be an efficient alternative to the mechanical 
ones. Development of the latter is restrained in particular by insufficient so far accuracy in 
damping identification. 
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FIG. 3: Imaginary and real parts of the eigenvalues of the gyroscopic system ([9]) as functions of 
the damping parameter 6i for fci = 1, = 0.3 and (a,b) k = 0, ^2 = —^i ("PT-symmetric case), 
(c,d) K = 0.1, 62 = -0.3. 

A gyroscopic system with indefinite damping. Taking into account commercial availability 
of q yrators — the non-reciprocal elements of LRC circuits that model gyroscopic effects 



45 



systems with the indefinite damping |l5 |. 



48| — it should be possible to extend the experiments described in j3] to the gyroscopic 



Consider a system with two degrees of freedom 

z + (D + 2^]J)z + (K + (^]J)2)z = 0, (9) 
where J is a matrix of gyroscopic forces with the entries ju = J22 = and J21 = — J12 = 1, 
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FIG. 4: (a) The domain of asymptotic stability and its boundary for the gyroscopic system ([9]) 
in the {k, X, y)-space when ki = 1 and O = 0.3. The vertical red interval of self-intersection 
corresponds to the domain of marginal stability of the P7~-symmetric gyroscopic system with the 
indefinite damping, (b) The top view of the stability boundary. 

is a gyroscopic parameter, and D and K are matrices of damping and potential forces. 
Eq. ([9]) describes stability of a particle in a rotating saddle trap and flexible shafts in the 
classical rotor dynamics and arises in the theories of helical quadrupole magnetic focussing 



systems of accelerator physics and light propagation in liquid crystals 22l. |49|-|52|. 



When D = diag {6i, —6i) and K = diag {ki, ki), the system iQ is invariant under trans- 



formations t ~t and Zi ^ Z2, i.e. it is PT-symmetric [t]. 

Assume D = diag 62) and K = diag {ki, ki + k). In Fig. [3] we plot the imaginary and 
real parts of the eigenvalues as functions of 61. When 62 = — 5i and k = 0, the spectrum 
is symmetric with respect to the imaginary axis of the complex plane and demonstrates a 
typical for the PT-symmetric system behavior, see Fig. [3t^a,b). Detuning the gain and loss 
as well as the potential, unfolds the EPs and creates an interval of the asymptotic stability 
that is smaller than the interval of the marginal stability, see Fig. [3]^c,d). 

With the parameters X = 61 + 62, and Y = 61 — 62, we plot the Routh-Hurwitz threshold 
for the asymptotic stabihty of system ([9]) in the (k, X, Y) space in Fig. HI Again, the surface 
is locally equivalent to the Pliicker conoid. PT-symmetric marginally stable systems live on 
the vertical interval of self-intersection terminated by two exceptional points. The Whitney 
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umbrella singularities at the EPs are responsible for the dissipation-induced enhancement 



of instability found in 



23| 



Example. Dissipatively enhanced modulational instability 

A monochromatic plane wave with a finite amplitude propagating in a dispersive medium 
can be disrupted into a train of short pulses when the amplitude exceeds some threshold. 
This process develops due to an unbounded increase in the percentage modulation of the 
wave, i.e. instability of the carrier wave with respect to modulations. This is a fundamental 



for modern fluid dynamics, nonlinear optics and plasma physics modulational instability |53 |. 



This instability, discovered by Bespalov and Talanov and Benjamin and Feir 



54 



55| . can 



trigger formation of the breather-type solitons from the Stokes waves in deep water. The 



breathers are associated with the rogue waves, recently detected in a water wave tank [56|. 

The modulational instability can be enhanced with additional dissipation [23] . Below we 
show that this effect is rooted in the mutual location of PT-symmetric gyroscopic systems 
with indefinite damping with respect to general dissipative ones. 

Without dissipation, a slowly varying in time envelope A of the rapidly oscillating carrier 
wave is often described by the nonlinear Schrodinger equation (NLS) 

iAt + aA^^ + -f\A\^A = 0, (10) 

where a and 7 are positive real numbers, i = a/— T, and the modulations are restricted to 



one space dimension x 



23 



53| . Eq. (|T0|) has a solution in the form of a monochromatic wave 
A = Aoe**^^-''^*, (11) 



where the frequency of the modulation, u, depends on the amplitude Aq = + iw^ and 
spacial wavenumber k as u = ak^ — 7||uo||^ with Uq = (m^,^^)- 

We linearize the NLS about the basic traveling wave solution f lTTj) in order to study 
stability of the modulation. Assuming periodic in x perturbations with the wavenumber a 
we substitute their Fourier expansions into the linearized problem. Then, the cr-dependent 
modes decouple into four-dimensional subspaces for each harmonic with the number n, so 



that for n = 1 we get [23 1 



Jv + 2akaJw — aa^v + 27Uouj^v = 0, 
Jw — 2akaJv — aa'^w + 27UoUq w = 0, (12) 



10 




b) 



c) 



0.5- 




T 1 1 1 1 1 r- 

1 




Modulational ^ 
instability 



0.4 



FIG. 5: (a) The boundary between the domains of asymptotic stability and modulational instability 
in the (o,c, ||uo||) space when u = 1, a = 1, 7 = 1, and k = 1. (b) In (c, ||uo||) plane the 
cross sections of the boundary at (green) a = and (red) a = 0.1. (c) In (a, c) plane the cross 
sections of the boundary at (green) ||uo|| = ||uo||i = (black) ||uo|| = ||uo||j — 0.05, and (red) 
||uo|| = ||uo||j + 0.05; blue lines is a linear approximation p8]) . 



where dot indicates time differentiation and the dyad uqu^ is a 2 x 2 symmetric matrix. 
Eq. f[T^ can be transformed to that of the indefinitely damped gyroscopic system with 
= acr^ — 7||uop, D = 27(uoUq J — JuqUq ), and K = {Aa'^k'^a'^ + 7^||uo||^)I, where I is 
a unit matrix, which is PT-symmetric because the eigenvalues ±27||uo||^ of the matrix D 



differ by sign only 22(|. This implies that the spectrum of the system (fT2|) is Hami 



i.e. symmetric with respect to both real and imaginary axis of the complex plane H 
with the eigenvalues 

A = ±i2aka ± ia\/ 2«7(||uo||^ — ||uoP), 



tonian, 



23 



where 



|Uo| 



aa 



27 



(13) 



(14) 



At small amplitudes of the modulation, the eigenvalues are pure imaginary. With the in- 
crease in the amplitude, the modes with the opposite Krein signature collide at the thresh- 



old ||uo| 



|Uo| 



23|. At ||uo|| > ||uo||t the double pure imaginary eigenvalue splits into 



complex-conjugate eigenvalues, one of which with positive real part, that corresponds to the 
modulational instability in the ideal (undamped) case 



23 



51. 



Introducing into Eq. f lTOl) the dispersive and nonlinear losses with the coefficients a and 



11 



c, respectively, we arrive at the dissipatively-perturbed NLS |23|-|25| 



iAt + (a - ia)A^^ + (7 + ^c)|ApA = 0, 



(15) 



which after hnearization and use of Fourier expansions yields the reduced system 23 1 



Jv + 2akaJw — aa v + 27UoUq v + 2kaaw + aa Jv + 2cUq vJuq = 0, 
Jw — 2«A;crJv — aa w + 27UoUq w — 2kaav + aa Jw + 2cUq wJuq = 0. (16) 

When a = and c = 0, Eqs. (fT6|) coincide with the ideal system (fT2l) . 

Writing the Routh-Hurwitz conditions for the characteristic polynomial of the system 
(fT6|) . we find an expression for the threshold of the modulational instability in the presence 
of dissipation 



2c {ca - 7a)||uo|| + {Aa ca{ca - 7a) - Aa'^k'^{Y + c ) + c a^{a^ + a 



2aa\aa\ac - 7a) + 2a^ca^ + Aak\'^a - ca))||uo|p + aV^(a^ - Ak^){a^ + a^) 



+ (17) 
= 0. 



The threshold equation (fTTI) yields a linear approximation to the stability boundary in the 
(a, c) plane of the coefficients of dispersive and nonlinear losses [23 1 



a 



-a ± 



A;(2||uo| 



|uo||iA/ll"o|li - ||uo| 



a + o(a). 



(18) 



When a <^ c, a simple approximation follows from Eq. (|T8|) to the amplitude at the threshold 
of the modulational instability in the presence of dissipation 

1 k^a^ c? 



|Uo| 



2||uo||f c2 



< lluol 



(19) 



Note that Eq. ( fT9ll is in the canonical for the Whitney umbrella form Z = X'^/Y'^ 41|. 

In Fig. [5t^a) the threshold f|T7|) is shown in the (a, c, ||uo||) space. At ||uo|| = ||uo||j and 
a = and c = it has the Whitney umbrella singularity at the exceptional point; along the 
interval ||uo|| < ||uo||j the system is PT-symmetric with pure imaginary spectrum. Below 
the surface (JT71) when a > and c > the dissipative system ( IT6l) with the broken VT- 
symmetry is asymptotically stable. In Fig. Mjo) the cross-sections of the stability boundary 
(ITTll are shown for a = (green line) and a = 0.1 (red line) in the (c, ||uo||) plane. The 
domain of modulational instability that was above the green line in Fig. [5](b) when a = 
expands considerably below this line (grey area) when the coefficient of dispersive losses 
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a 7^ {enhancement of the modulational instability with dissipation 23|]). Fig.[5]^c) shows the 
cross-sections of the surface f|T71) in the (a, c) plane for ||uo|| = ||uo||i (green hne) and when 
||uo|| is shghtly above (red hne) or below (black line) the amplitude at the threshold of the 
modulational instability in the undamped case. The cross-sections are typical for the surface 
with the Whitney umbrella singularity 32|]. In particular, they justify the approximation 
( !T8|) (blue lines) to the stability boundary that yields the canonical equation for the Whitney 
umbrella f lT9|) . 

Summary. A direct link is established between the PT-symmetry and dissipation-induced 
instabilities: The systems with the exact PT-symmetry are identified with the singularities 
on the threshold of asymptotic stability of the indefinitely damped ones. This finding opens 
a new perspective for PT-symmetric LRC circuit experiments that could test the both 
fundamental physical concepts, which is so far unavailable in the mechanical experiments. 
As an example, the enhancement of the modulational instability with dissipation is connected 
to the existence of the Whitney umbrella singularity on the instability threshold. 



[1] C. M. Bender and S. Boettcher, Phys. Rev. Lett. 80, 5243 (1998). 
[2] C. M. Bender, S. Boettcher, P. N. Meisinger, J. Math. Phys. 40, 2201 (1999). 
[3] A. Mostafazadeh and A. Batal, J. Phys. A-Math. Gen. 37, 11645 (2004). 
[4] C. M. Bender, Rep. Prog. Phys. 70, 947 (2007). 

[5] A. Guo, G. J. Salamo, D. Duchesne, R. Morandotti, M. Volatier-Ravat, V. Aimez, G. A. 

Siviloglou, D. N. Christodouhdes, Phys. Rev. Lett. 103, 093902 (2009). 
[6] C. E. Ruter, K. G. Makris, R. El-Ganainy, D. N. Christodouhdes, M. Segev and D. Kip, Nat. 

Phys. 6 192 (2010). 

[7] J. Schindler, A. Li, M. C. Zheng, F. M. Ellis and T. Kottos, Phys. Rev. A, 84 (2011), 
040101(R). 

[8] Q. Wang, Czech. J. Phys., 54, 143 (2004). 

[9] C. M. Bender and P. D. Mannheim, Phys. Lett. A, 374, 1616 (2010). 
[10] P. Freitas, M. Grinfeld and P. A. Knight, Adv. Math. Sci. Appl. 17, 435 (1997). 
[11] P. Freitas, ZAMP 50, 64 (1999). 

[12] G. Chen, S. A. Fuhing, F. J. Narkowich, and S. Sun, SIAM J. Appk Math. 51, 266 (1991). 



13 



[13] P. Freitas, J. DifF. Eqns. 127, 320 (1996). 

[14] R. Joly, Asymp. Anal. 53, 237 (2007). 

[15] W. Kliem and C. Pommer, ZAMP 60, 785 (2009). 

[16] O. N. Kirillov, Phys. Lett. A 373, 940 (2009). 

[17] T. Damm and J. Homeyer, On indefinite damping and gyroscopic stabilization, in: Prepr. 

18th IFAC World Congr., 7589 (2011). 
[18] R. T. Spurr, Wear 4, 150 (1961). 
[19] A. Akay, J. Acoust. Soc. Amer. Ill, 1525 (2002). 
[20] O. N. Kirillov, Proc. R. Soc. A. 464, 2321 (2008). 
[21] O. N. Kirillov, Proc. R. Soc. A. 465, 2703 (2009). 
[22] O. N. Kirillov, Phys. Lett. A 375, 1653 (2011). 
[23] T. J. Bridges and F. Bias, Phys. Fluids. 19, 104104 (2007). 
[24] B. A. Malomed and H. G. Winful, Phys. Rev. E 53, 5365 (1996). 
[25] J. Atai and B. A. Malomed, Phys. Rev. E 54, 4371 (1996). 

[26] A. E. Miroshnichenko, B. A. Malomed, and Y. S. Kivshar, Phys. Rev. A 84, 012123 (2011). 
[27] R. Driben and B. A. Malomed, Europhys. Lett. 96, 51001 (2011). 
[28] F. K. AbduUaev, Y. V. Kartashov, V. V. Konotop, D. A. Zezyulin 83, 041805 (2011). 
[29] O. N. Kirillov, Int. J. Non-Lin. Mcch. 42, 71 (2007). 

[30] A. M. Bloch, P. S. Krishnaprasad, J. E. Marsden, and T. S. Ratiu, Annales Inst. Henri 

Poincare, 11(1), 37 (1994). 
[31] O. N. Kirillov, Dokl. Phys., 49(4), 239 (2004). 
[32] O. N. Kirillov and F. Verhulst, ZAMM 90, 462 (2010). 

[33] B. Dietz, H. L. Harney, O. N. Kirillov, M. Miski-Oglu, A. Richter, F. Schafer, Phys. Rev. Lett. 

106 150403 (2011). 
[34] N. P. Hoffmann and L. Gaul, ZAMP 83, 524 (2003). 
[35] R. S. MacKay, Phys. Lett. A 155, 266 (1991). 
[36] R. Romea, J. Atm. Sci. 34, 1689 (1977). 

[37] R. Krechetnikov and J. E. Marsden, Arch. Rat. Mech. Anal. 194, 611 (2009). 

[38] G. E. Swaters, J. Phys. Oceanogr. 40, 830 (2010). 

[39] I. Hoveijn and M. Ruijgrok, Z. angew. Math. Phys. 46, 384 (1995). 

[40] I. Hoveijn and O. N. Kirillov, J. DifF. Eqns. 248, 2585 (2010). 

14 



[41] M. Berger and B. Gostiaux, Differential Geometry: Manifolds, Curves and Surfaces, Grad. 

Texts in Math., vol. 115, Springer, Berlin (1988). 

[42] O. Bottema, Indag. Math. 18, 403 (1956). 

[43] S. A. Van Gils, M. Krupa, W. F. Langford, Nonlinearity 3, 825 (1990). 

[44] W. F. Langford, Hopf meets Hamilton under Whitney's umbrella, in: lUTAM Symp. Nonlin. 

Stoch. Dyn., Solid Mech. Appl. 110, Kluwer, Dordrecht, 157 (2003). 

[45] B. D. H. Tellegen, Philips Res. Kept. 3, 81 (1948). 

[46] A. Fabre, Electron. Lett. 28(3), 263 (1992). 

[47] A. Figotin and I. Vitebsky, SIAM J. Appl. Math. 61, 2008 (2001). 

[48] R. J. Potton, Rep. Prog. Phys. 67, 717 (2004). 

[49] L. E. J. Brouwer, N. Arch. v. Wisk. 2, 407 (1918). 

[50] O. Bottema, ZAMP 27, 663 (1976). 

[51] C. W. Roberson, A. MondelH and D. Chernin, Phys. Rev. Lett. 50, 507 (1983). 

[52] V. E. Shapiro, Phys. Lett. A 290, 288 (2001). 

[53] V. E. Zakharov and L. A. Ostrovsky, Phys. D 238, 540 (2009). 

[54] V. Bespalov and V. Talanov, JETP Lett.-USSR 3, 307 (1966). 

[55] T. B. Benjamin and J. E. Feir, J. Fluid Mech. 27, 417 (1967). 

[56] A. Chabchoub, N. P. Hoffmann and N. Akhmcdicv, Phys. Rev. Lett. 106, 204502 (2011). 



15 



